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ABSTRACT 

To study the current structure of the California Current as it manifests itself 
on the continental shelf a subsurface mooring, PI, was anchored 5 km west of 
Point Sur at 36° 17' N, 121° 59' W from 28 February through 1 1 May 1990. The 
PI mooring, placed on the 84 m isobath, consisted of a Self-Contained Acoustic 
Doppler Current Profiler (SC-ADCP) housed in a syntactic foam sphere and 
secured to an anchor. The mooring geometry placed the transducer heads at 80 
m depth. The instrument operated at 307 kHz with each beam inclined 30° from 
the vertical. An ensemble average of 170 one-second pings were recorded every 
15 minutes for 70 days. Good velocity data to within 15 m of the surface were 
obtained at 4 m intervals throughout the record. 

Low-pass filtered data indicate high coherence between the observed 
currents and the local wind forcing measured at the Monterey Bay NOAA buoy 
located 53 km to the northwest of PI. During equatorward wind events the 
current response was vertically sheared and time-lagged at depth with near 
surface magnitudes up to 35 cm s'^. During wind relaxations (periods of no 
wind) and poleward winds a pressure induced barotropic poleward flow, on the 
order of 10 cm s'^, developed throughout the water column. These poleward 
currents became more frequent and better developed later in the record. Analysis 
of adjusted sea level (ASL) differences between stations located to the north and 
south of mooring PI suggest an increasing poleward pressure gradient during the 
span of the record. Alongshore differences in the alongshore wind stress 
component contributed to the increasing difference in ASL. 

A simple wind forced theoretical model adapted from Ekman (1905) and 
Csanady (1982), was used to describe the flow in the coastal regime. The model 
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consists of surface and bottom boundary layers superimposed on a mostly 
geostrophic interior flow driven by set-up and set-down of the sea surface near 
the coast. The magnitude, direction, and sense of rotation of the current vectors 
with depth were reproduced well by the model when it was forced by the actual 
magnitude of the surface wind stress and a realistic vertical eddy coefficient (Az). 
Differences between the model and the observations can be accounted for by 
other elements of the coastal circulation not included in the model. 
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II. INTRODUCTION 



A. THE CALIFORNIA CURRENT SYSTEM 

The California Current System (CCS), an eastern boundary current system, is 
characterized by a flow that is highly variable both temporally and spatially. The 
CCS is composed of the following components: the California Current (CC), the 
California Undercurrent (CUC), the Davidson Current (DC) and the Southern 
California Counter Current (SCCC). The CC, comprising the eastern limb of the 
north Pacific anticyclonic gyre, flows along the west coast of the United States 
and is characterized by predominantly equatorward mean flow with velocities on 
the order of 1-10 cm s‘^. The majority of the equatorward transport occurs within 
1000 km of the coast with the current reaching from the surface down to 300 m 
depth (Hickey, 1979). Flow over the continental shelf tends to be highly variable 
with equatorward and poleward surface currents superimposed on the mean flow 
(Chelton, 1984). At depths greater than 100 m, on the continental slope, the CUC 
is present throughout the year and consists of a poleward component with 
seasonally averaged velocities on the order of 10 cm s'^ (Chelton, 1984) and 
maxima up to 40 cm s'^ (Tisch, 1990). Only present in the fall and winter months, 
the DC is a poleward surface current north of Point Conception occurring within 
80 km of the coast with a magnitude on the order of 10 cm s' ^ (McCreary et 
al., 1987). The companion to the DC is the SCCC which is also poleward and 
encompasses the area south of Point Conception (Hickey, 1979). The CCS is 
generated by large scale atmospheric forcing between the north Pacific high and 
the continental thermal low over the southwest U. S. The seasonal variation of 
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the winds is evident in Figure 1. The high pressure migrates from its most 
southerly position of 30° N, 135° W in the winter months to its most northerly 
position 38° N, 145° W in the summer. There is a corresponding building of the 
high in the summer months (1020 mb in November to 1025 mb in August). The 
low pressure over the southwest portion of the U. S. does not tend to migrate 
seasonally, however, there is a deepening in the summer months due to warming 
of the land (Huyer, 1983). The pressure gradient associated with these features 
then will also exhibit seasonal variation. 




Figure 1. Geographical position of seasonal atmospheric high (H) and 
low (L) pressure systems over the eastern North Pacific and the west coast of 
North America. Isoheight pressure levels are contoured with the direction of 
flow indicated by arrows. Pressure is measured in mb (from Reid et al., 1958). 
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As the pressure gradient manifests itself by the wind field, in the summer months 
south of 40° N winds are generally steady, equatorw'ard, and highly coherent in 

the alongshore direction (Chelton, 1984). 

Steady equatorward winds produce offshore flow via the Ekman transport 
which generates dynamic sea level height differences perpendicular to the coast. 
This results in a geostrophically balanced surface flow parallel to the coast 
located typically within 5 - 25 km from the coast which typically encompasses 
the internal Rossby radius of deformation (Hickey, 1979). Since variations in the 
current regime respond to local wind forcing, changes can occur rapidly (on the 
order of days). The magnitude of the wind forcing is strongest approximately 
200 km off shore. As a result, positive wind stress curl (curl x) is present in the 
near shore region (Nelson, 1977), and may be important in driving poleward flow 
(McCreary et al., 1987). 

B. COASTAL UPWELLING 

Coastal upwelling occurs seasonally along the entire west coast of the U. S. 
and is the result of equatorward winds generating offshore Ekman transport 
away from the coast (Huyer, 1983). Most of the knowledge regarding the 
structure and dynamics of the active upwelling zone is limited to a coastal strip a 
few tens of kilometers wide (Huyer, 1983). To help illustrate the dynamics of 
coastal upwelling circulation on the continental shelf a conceptual two 
dimensional model (no alongshore variation) is presented (Figure 2). 
Equatorward winds drive surface water off-shore in the Ekman layer. Below this 
layer, a compensating return drift is generated to conserve volume. This return 
flow will then surface in a band near shore. Near the surface, the surfacing 
thermocline produces a baroclinic jet, whose vertical shear is geostrophically 
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balanced by the across-shore density gradients induced by upwelling (Kosro, 
1987). At greater depths there is a reversal of the along-shore currents, a 
poleward undercurrent, which flows counter to the wind (Kosro, 1987) 



Figure 2. Schematic of two dimensional upwelling. (a) This figure 
represents the cross-shore (u) component of flow. Current speed and direction 
are indicated by vectors. The surface and bottom Ekman layers are indicated by 
dashed lines, (b) This figure represents the alongshore (v) component of flow. 
Velocity contours are indicated with equatorward flow occurring in the surface 
layers and a poleward undercurrent present near the sea floor (from Huyer, 1983). 



In contrast to the two dimensional schematic, which is useful in describing 
the upwelling concept, in many applicable flow regimes it appears that coastal 
upwelling is always a three dimensional process (Brink, 1983), imposed by basin 
geometry or by external forcing imposed on a large scale. This results in certain 
areas ("upwelling centers") near coastal promontories where upwelling appears to 
be more vigorous than others. Based on satellite imagery (Tisch, 1990), Point Sur 
may act as an upwelling center. 




(a) 



(b) 
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C. PAST STUDIES 

Extensive study in the large scale circulation of the CCS was initiated in 
1949 by the California Cooperative Fisheries Investigation (CalCOFI) to examine 
the ecosystem of the southern portion of the CCS (Chelton et al., 1988). An 
important aspect of this CalCOFI study was the maintenance of a geographically 
fixed sampling grid to take hydrographic surveys. This sampling grid consisted of 
parallel lines oriented perpendicular to the coast and covered a region from San 
Francisco to Southern Baja. With a typical cross-shelf grid spacing of 74 km, 
large-scale low-frequency variability in the current structure is well addressed 
(Chelton et al., 1988). However, detailed temporal and spatial analysis of smaller 
scale high-frequency events occurring along the coastal regions cannot be 
obtained from this study. 

1. Coastal Oceans Dynamics Experiment (CODE) 

The Coastal Oceans Dynamics Experiment (CODE) was conducted 
along the continental shelf from Point Reyes to Point Arena in the spring and 
summer of 1981 (CODE 1), and 1982 (CODE 2). The goal of CODE was to obtain 
high quality data sets needed to dynamically describe the response of the 
continental shelf waters to wind forcing in the temporal band covering the 2-20 
day variability (Beardsley and Lentz, 1987). 

Atmospheric forcing during CODE 1 and 2 was characterized by periods 
of strong (7-15 m s'^) upwelling favorable alongshelf winds lasting up to 30 d, 
interrupted by shorter periods of much weaker winds directed either poleward or 
equatorward. The periods of weak or reversed winds typically lasted several 
days and were primarily associated with coastally trapped perturbations of the 



5 



marine layer along the central California coast and not weakening of the North 
Pacific high. (Beardsley et al., 1987) 

Maximum low’ frequency response in the water column to equatorward 
wind forcing was found at midshelf and near surface. The mean current during 
the summer upwelling season was strongly sheared in both vertical and 
horizontal cross shelf directions. The flow over the outer shelf was surface- 
intensified and directed off-shelf and equatorward. The flow over most of the 
water column over the inner shelf was directed poleward. In the absence of local 
wind forcing, a barotropic jet formed with the maximum amplitude forming over 
the inner continental shelf. The effect then, of strong equatorward wind forcing 
would be to counter the poleward flow so that currents on the shelf are 
everywhere equatorward. (Winant et al., 1987). 

The coastal current was measured by a shipboard Doppler acoustic log 
during CODE (Kosro, 1987) and provided detailed information on the spatial 
structure of the current field. The nearshore current field was found to be 
anisotropic, varying more rapidly across-shore than alongshore. A surface- 
intensified equatorward jet was found in the average alongshore currents over 
much of the shelf. At depths greater than 80 m the average alongshore flow 
reversed producing a poleward undercurrent, with strongest amplitudes near the 
shelf break. Mean current fluctuations were strongly polarized along isobaths 
near the coast. (Kosro, 1987) 

Two other CODE papers indicate the importance of remote forcing from 
the south along the California coast (Denbo and Allen, 1987; Davis and Bogden, 
1989). The basic tenet in either case is that the local current dynamics are 
consistent with wind-forced, long, coastal trapped waves which are generated to 
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the south and propagate northward in the coastal waveguide. Denbo and Allen 
(1987) find that the short period (< 4 - 7 d) are more locally forced, while the long 
periods (> 4 - 7 d) are best correlated with the wind stress 170 - 300 km to the 
south. After removing the newly-discovered "Kelvinlike" wave from the sea level 
and bottom pressure observations, Davis and Bogden (1989) found that the 
maximum coherence between observed currents on the CODE site and coastal 
wind stress was with winds measured 500 km to the south. 

2. The Central California Coastal Circulation Study (CCCS) 

To gain more understanding of the circulation on the continental shelf 
between Point Conception (PC) and San Francisco (SF), the Central California 
Coastal Circulation Study (CCCS) was conducted February 1984 through July 
1985. Current meter moorings and high resolution hydrographic surveys with a 
cross-shore spacing of 20 km were conducted inshore of the continental shelf 
break to study the seasonal variability and spatial coherence of the flow over the 
continental shelf (Chelton et al., 1988). 

From February 1984 through July 1984 there was strong mean 
poleward flow on the order of 20 cm s'^ at the 70 m level along the shelf between 
Point Sur (PS) and PC in opposition to equatorward wind stress with the mean 
flow north of PS which was equatorward. Fluctuations in alongshore currents 
over the shelf were highly correlated with the wind everywhere except PC. 
During periods of weak equatorward winds the poleward continental shelf flow 
south of PS spread further offshore to the continental slope. Similar nearshore 
flow had previously been observed off Point Sur throughout 1979 and 1980 
(Wickham et al., 1987). The results of CCCS suggest that anomalous weak winds 
along the central California coast in June and July 1981 and 1984 may have been 
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responsible for the intensification of poleward flow on the continental shelf 
(Chelton et al., 1988). 

D. PURPOSE OF STUDY 

The Naval Postgraduate School (NPS) Department of Oceanography 
maintains current meter moorings at 800 m and 1800 m along 36° 20' N to 
measure the energetics of the CC and the CUC as they flow along the central 
California continental slope (Figure 3). The reader is referred to Sielbeck, (1991) 
for a complete description of these moorings and a discussion of the current meter 
accuracies and deployment procedures. To gain a better understanding of how 
the flow is coupled to the continental shelf a Self-Contained Acoustic Doppler 
Current Profiler (SC-ADCP) was deployed on the 84 m isobath from 28 February 
through 11 May 1990. The original intent was to place this mooring in line with 
the others, however, due to fishing activity the mooring was placed at 36° 17' N, 
121° 59' W approximately 5 km west of Point Sur. 

The SC-ADCP over the past several years has been proven as a reliable 
means to obtain vertical profiles of horizontal ocean currents (Schott and Johns, 
1987; Schott, 1986; Johns, 1988; Pettigrew, 1988). The SC-ADCP allows the user 
the flexibility to place the instrument in the relatively secure surroundings of the 
sea floor and still observe current in the upper water column. The purpose of this 
study is to observe the current structure off Point Sur, California using an SC- 
ADCP and relate these currents to the local and remote forcing. Data acquisition 
techniques are described in Chapter II. Data processing and supplemental data 
sets are discussed in Chapters III and IV. The current data set observations are 
presented in Chapter V and a theoretical coastal current model, forced by the 
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observed winds, is described in Chapter VI. Discussion, conclusions and 
recommendations conclude this paper. 



- 122.50 - 122.25 - 122,00 - 121.75 - 121.50 




Longitude 



Figure 3. Chart of the California coast off Point Sur locating the mooring 
positions of moorings PI, P2, P3: P3 is located on the 1800 m isobath, P2 is 
located on the 800 m isobath, and PI is located on the 84 m isobath. 
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II. DATA COLLECTION 



A. PRINCIPALS OF OPERATION 

L The SC-ADCP and mooring geometry 

The SC-ADCP (Figure 4a), manufactured by RD Instruments Inc., is a 
modular instrument consisting of a water-tight transducer assembly which 
operates at an acoustic frequency of 307 KHz, and a pressure case that houses a 
magnetic tape recorder, battery packs, pitch and roll sensors, fluxgate magnetic 
compass, and a high resolution temperature sensor. The SC-ADCP was deployed 
from 28 Feb - 1 1 May 1990 at site PI. The mooring was in 84 m of water on a 
relatively flat section of the sea floor. A rocky shoal area 55 m deep was located 
about 500 m to the southeast (Figure 3). 

The SC-ADCP was secured in a 49 inch syntactic foam sphere that 
provided approximately 1000 pounds of net buoyancy with the SC-ADCP in 
place (Figure 4b). The sphere was attached to an anchor via a chain assembly 
and dual acoustic releases. Because of the large net buoyancy force, the mooring 
acted essentially as a rigid body experiencing very little tilt (< 1°) or rotation 
throughout the deployment. The length of the mooring, approximately 4 m, 
placed the transducer heads at the 80 m depth. 

2. Velocity calculation 

The SC-ADCP operates on the Doppler principal. Acoustic pulses are 
transmitted along the four vertically inclined beams and the backscattered 
frequency shifts from scatterers within the water column are measured at the 
transducer head (Johns, 1988). Since only radial motion relative to the transducer 
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Figure 4. SC-ADCP configuration and PI mooring geometry: (a) The SC- 
ADCP is a modular instrument consisting of a transducer assembly with four 
mutually orthogonal transducers at an angle of 30® from vertical and a pressure 
case that houses the electronics, (b) The SC-ADCP is housed in a syntactic foam 
sphere that provides 1000 pounds of net buoyancy. 
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head will cause a Doppler shift, only the velocity component in the direction of 
the line between the source and the receiver will be observed. The Doppler shifts 
of the beams are given by 



r 2Vj{ft sin 6) 
iDi = ;; 



i=l,2.3.4 



( 1 ) 



where fpj is the frequency shift along each beam, Vi is the relative velocity 
between source and receiver, fj is the transducer transmit frequency, c is the speed 
of sound at the transducer head, and 6 is the angle of the transducer head from 
the vertical. 

Equation (1) is solved to find the radial velocity component (Vi) as a 
function of the observed Doppler shift (fpj) and further post processing is 

required to convert the radial velocity component into earth-based velocity 
components. Programming options are available that cause the SC-ADCP to 
automatically transform the radial velocities into earth coordinates for each ping, 
to allow in-situ vector averaging. Utilizing this option, and the known beam 
geometry, radial velocities from three beams are converted into east (u), north (v), 
and vertical (w) velocity components. The fourth beam is redundant and is used 
for error checking (RD Instruments, 1989). 

B. INSTRUMENT CONFIGURATION 

The SC-ADCP was configured to record an ensemble average of 170 acoustic 
pulses, or pings, every 15 minutes with a ping duration of one second. The SC- 
ADCP segmented the velocity profile into uniform depths called depth cells or 
bins by range-gating the return from each ping. Twenty 4-m bins were recorded 
for each of the 7025 usable ensembles (record length of 70 days). Adjacent cells 
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were not statistically independent, because of the overlap from the sliding digital 
filter which caused a correlation between adjacent cells of 0.15 (RD Instruments, 
1989). In order to examine statistically independent current data, only every 
other bin was used for further processing. This in effect produced a vertical 
profile consisting of 10 independent current measurements centered at 8 m apart. 
The stored data consisting of header and current information were converted into 
ASCII format utilizing software provided by RD Instruments. The current data 
were then transferred to the NFS Vax 11/780 computer where the count 
formatted data were converted into engineering units. All subsequent analysis 
and graphical display was accomplished using the Vax computer. 

C. ERROR ANALYSIS 

1. Bias and random error 

The nature of Doppler velocity measurement will induce two kinds of 
errors that contribute to velocity uncertainty. These errors are bias and random 
error. Bias is typically on the order of 0.5 - 1.0 cm s' The bias of the instrument 
depends on a variety of factors including temperature, mean current speed, 
signal/noise ratio, beam geometry errors etc. It is not yet possible to exactly 
measure SC-ADCP bias and to remove it in post-processing. 

Single ping SC-ADCP operation will also introduce random error. Since 
random error is uncorrelated for each ping, averaging multiple pings will reduce 
the standard deviation of the error velocity. Because the external bias of the SC- 
ADCP is on the order of 0.5 - 1.0 cm s'^ it is not beneficial to average the ping 
data below this value. The number of pings per ensemble that must be averaged 
to reduce the standard deviation of the random error to a specified level is given 
by 
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( 2 ) 



p^( 1 .^x1qM 2 

where P is the number of pings per ensemble, F is the SC-ADCP transmission 
frequency, SD is the desired standard deviation in cm s"^, L is the depth cell 
length, and B is the beam angle coefficient (1.0 for a 30° beam). Equation (2) 
yields a 170 ping ensemble necessary to resolve the standard deviation of random 
error velocity to 1.0 cm s'^ which is equal to the bias of the instrument. (RD 
Instruments, 1989) 

2. Sound speed difference at transducer 

The SC-ADCP assumes a fixed sound speed of 1536 m s'^ over the 
entire water column. Since the recorded velocity values are scaled frequency 
shifts in Hertz, horizontal velocity errors will be introduced if there are variations 
in sound speed. To correct for variations in sound speed a constant 
multiplicative factor of Co / 1536 m s‘l is used where Co is the speed of sound at 
the transducer head. This correction is only required at the transducer head 
because variations in sound speed over the water column are exactly 
compensated for by ray bending (Johns, 1988). Snell's Law stipulates that 
horizontal wavenumber must be conserved when sound passes through 
horizontal surfaces. Since the frequency is constant, sound speed variation, as a 
function of depth will not affect the horizontal component of the sound velocity 
(RD Instruments, 1989). 

3. Conversion to earth coordinates 

The velocity data retrieved from the SC-ADCP deployed off Point Sur 
were stored on a magnetic tape in earth-based reference coordinates. The internal 
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processing required to convert the radial-coordinate beam data into earth- 
coordinate based system is described with emphasis placed on sources of 
potential error. 

When the SC-ADCP is mounted on a mooring free to change orientation 
the automatic conversion of the radial velocity components to earth-based system 
is required. Recording the data in earth-coordinates reduces post-processing, 
although there may be a slight trade-off in data quality. The main loss can be 
attributed to the depth cell mapping feature of the instrument. The SC-ADCP 
automatically calculates velocity at particular depth cells that are at the same 
depth (RD Instruments, 1989). If excessive tilt is present, the SC-ADCP 
processing software will try to adjust for the deviation, from the normal 
orientation of straight up, by reassigning bins from each beam to different depth 
levels so that the velocity is always calculated from bins at about the same depth 
(Johns, 1988). This may cause some bins to be skipped. The observed tilt at 
mooring PI (Figure 5), however, indicates that the current induced mooring tilt 
was negligible and that errors introduced would be insignificant. 

D. DATA VELOCITY CHECKS 
1. Echo Amplitude 

Echo amplitude is a measure of the energy intensity of the backscattered 
echo for each depth cell. Useful in determining signal strength, it is also used as a 
check for data quality. Since echo amplitude is a function of the distance from 
the transducer and the number of scatterers in the water column, echo amplitude 
will normally decrease until the theoretical range of the instrument is reached. For 
the 307 KHz SC-ADCP this range is 120 m and the surface should be well within 
range of the instrument (RD Instruments, 1989). 



15 





1 


iiiiiiiiiiiiiiiijiiiiiiiiiiiiiiiiiiiijijiiiiiiijjiiiiiiiiiiiiiiiiiiiiii 




0 I 

1 1 

^ 90. - 


1 


iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii 


1 1 1 1 1 1 1 


/ U. 


1 


IlllllilllllllllillllllllllllillilllllJliJlililllllllJillilllllllllllll 




PITCH 

1 1 1 

r“ P 

k> o bo 

o o o 

1 1 1 1 1 1 






— 


iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiimmiimmmimimmmimmiimi 



28 7 14 21 28 7 14 21 28 7 

MAR APR MAY 



0.40 

0.20 

0.00 



o 

q: 



Figure 5. Time series of heading, pitch, and roll from the SC-ADCP: 
Heading is degrees magnetic, roll is in degrees and represents deviation in the u 
plane of motion, and pitch is in degrees and represents deviation in the v plane of 
motion. The constant offset of about 1® in pitch is due to the way the SC-ADCP 
was mounted in the buoy and was not induced by the current. 



The values of echo amplitude for all four beams at each depth level, 
averaged over the entire deployment (Figure 6) show that the amplitude of the 
signal decreases as a function of range from the transducer heads. However, in 
the near surface bins, there is a substantial increase in amplitude. The large 
amplitudes can be accounted for by the characteristics of the acoustic pulses. 
Each pulse consists of a main lobe where most of the power is directed and 
sidelobes which contain the remaining power. The sidelobes are designed to be 
at least 40 dB below the main lobe and are found about 40® offset from the main 
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average values of echo amplitude for all four beams averaged over the 
deployment with standard deviation superimposed over each value. Percentage 
of good pings are represented by percent bars and are averaged for all four 
beams (from RD Instruments, 1989). 
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lobe. The echo from a hard surface such as the sea surface is much stronger than 
the echo from the scatterers in the water. With each transducer oriented at 30° 
from the vertical, the echo from the side lobe facing the surface will return to the 
SC-ADCP at the same time as the echo from the main lobe at 85 percent of the 
distance to the surface (Figure 7), which overwhelms the side-lobe suppression 
and regulates the quality of the data. The equation that governs this is given by: 

Rmax = D cos 0 

where Rmax is the maximum range for acceptable data, D is the distance from the 
SC-ADCP to the surface, and 0 is the angle of the acoustic beam relative to 
vertical. For an SC-ADCP at 80 m depth with a 30° beam angle, the value of 
Rmax = 69 m and data within 1 1 m of the surface will be contaminated and should 
be disregarded. This is indicated by the increase in echo amplitude in the upper 
three bins (Figure 6). (RD Instruments, 1989) 

The low echo amplitude associated with bin 1, or the deepest level, can 
be attributed to a hardware problem in the receive circuitry of the SC-ADCP (RD 
Instruments, Personal Communication). Based on analysis covered in later 
sections, this data was determined to be usable. 

2. Percent Good 

The percent good is a measure of the percentage of pings per ensemble 
that are greater than the signal-to-noise threshold programmed into the SC- 
ADCP. The threshold was set to 6 dB above background noise (RD Instruments, 
1989). Each percent bar (0-99) (Figure 6) represents the percentage of good 
pings averaged over all four beams during the deployment of the instrument. 
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Values indicate high confidence in the data collected encompassing the water 
column. There is however, a difference of about seven percent near the surface 
probably due to side-lobe interference. 
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Figure 7. SC-ADCP transducer beam characteristics: The design 

specification for the side lobe is 40 dB below the main lobe and offset by 40®. 
(from RD Instruments, 1989) 
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III. SUPPLEMENTAL DATA SETS 



Wind and atmospheric pressure data from National Oceanographic and 
Atmospheric Administration (NOAA) weather buoys stationed off the California 
coast and sea level data from tide stations operated by the National Ocean 
Service (NOS) were utilized in this study. The geographical locations of these 
instruments are summarized in Figure 8. 

A. WIND STRESS DATA 

Wind data was available from the Monterey Bay buoy and the Cape San 
Martin buoy located to the north and south of site PI respectively. Strong 
correlation between the data sets suggested local wind forcing at PI could be 
taken from either source. The local stress was derived from wind data at the 
Monterey Bay buoy following the method of Large and Pond (1981). The reader 
is referred to Tisch (1990) for gap filling techniques and wind level correction 
schemes employed. Once all data gaps were filled and wind speed corrected to 
the 10 m level, surface stress components were calculated according to 

tT = Pa CdM u 

V = Pa CdM V 

where Xy and x* are the northward and eastward components of wind stress, pa 
is the atmospheric density, Cd is the drag coefficient, and M is the magnitude of 
the velocity with the u - component east and the v - component north. 
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Figure 8. Chart of California coast indicating geographic positions of 
NOAA buoys and tide stations. 
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B. HYDROGRAPHIC DATA 

TTie determination of the water characteristics, particularly the stratification, 
was essential to accurately model the flow field. Several hydrographic surveys 
were conducted in the vicinity of mooring PI in support of related activities of 
the NPS Oceanography Department, Data were collected from the surface to the 
near bottom using a Neil Brown Instruments System Mark III CTD. The reader is 
referred to Tisch (1990) for a detailed description of CTD data acquisition 
techniques and procedures. 

C. SYNTHETIC SUBSURFACE PRESSURE 

Synthetic subsurface pressure (SSP), also known as adjusted sea level (ASL), 
was calculated as the algebraic sum of the atmospheric pressure and the local sea 
level converted to decibars (Brown et al., 1985). This procedure partially corrects 
for the inverse barometer effect so that the remaining pressure fluctuations are 
primarily due to oceanic phenomenon. The calculation was made at two coastal 
sites, Monterey Bay and Port San Luis. Atmospheric pressure used for Monterey 
Bay was based on the average of Monterey Bay and Cape San Martin, and for 
Port San Luis, the average of Santa Maria and Port San Luis. The reader is 
referred to Tisch, (1990) for gap filling procedures and data conversion 
techniques. To obtain the alongshore difference in SSP the two-year mean value 
was removed for each station and the Monterey Bay data was subtracted from 
the Port San Luis data. The sign convention is such that a (-t-) SSP difference 
implies a poleward pressure gradient and a (-) SSP difference, an equatorward 
pressure gradient. Since the two stations were not geophysically leveled, the true 
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mean difference in SSP height is unknown, however, trends and fluctuations in 
ASSP should be well represented by the technique. 
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IV. DATA PROCESSING 



A. DATA FILTERING 

After all data was downloaded into the Vax computer each data set was 
filtered using methods originally developed by Oregon State University (Denbo 
et al., 1984). Current data at 15 minute intervals was first low passed (LP) filtered 
using a Cosine-Lanczos filter. This filter utilized a one half power period of 2.9 
hours (or 8.4 cpd) which allowed any signal with a period exceeding 3.6 hours to 
pass unattenuated. A centered 25 point segment of the record was used to 
produce each LP data point. After LP filtering each record was interpolated to 
transform 15 minute data records to 60 minute intervals with each data point 
falling on a even hour. To eliminate tidal signals, the current data was then low 
pass filtered (LLP) using a Cosine-Lanczos filter with a constant 121 point 
segment and had a half power period of 46.59 hours. The data were subsampled 
to transform the 60 minute data to 6 hour intervals with each data point falling on 
an even 6 hour mark. (Sielbeck, 1991) 

Since the wind data were recorded at 60 minute intervals, only the LLP 
filtering was accomplished. The data were then subsamplcd to six hour intervals 
as with the current data. A sample from the SC-ADCP data set (Figure 9) shows 
the result of filtering. 
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Figure 9. Example of filtering technique: Raw data for each component is 
passed through a low pass (LP) filter and plotted at hourly marks. Hourly data is 
then passed trough a low pass (LLP) filter and plotted at six-hourly marks. 
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B. SPECTRAL ANALYSIS 



To quantify the graphical presentations in Chapter IV, spectral analysis was 
performed using the spectral processing package PROSPECT developed by the 
Woods Hole Oceanographic Institution, This software package was written in 
Fortran 77 and designed to run on the VAX 11/780 computer. It performs 
spectral analysis of equally-spaced time series. 

The following discussion details the basic operations required to estimate 
coherence and phase functions. This discussion utilizes theory following Bendat 
and Piersol (1986) and operations performed by PROSPECT. The current and 
wind stress data sets consist of 283 equally spaced points of six-hourly LLP data. 
The data sets were divided into three pieces of 92 points each. To suppress 
sidelobe leakage, a Hanning taper was applied to each block of data. The Fast 
Fourier Transform (FFT) was computed for each block and in order to reduce the 
variance increase caused by tapering, the FFT was computed for a 50 % 
overlapped record. 

The auto spectra density estimate is given as 

Gxx(fk) = Z I Xi(fk)P k = 0. 1, ...| 

where n is the number of blocks of data, N is the number of data values, and At is 
the time step. The cross spectral density estimate for each pair of blocks is 
computed as 



Gxy(fk) = Cxy(fk) - iQxy(fk) 
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where Cxy(fk) is the cospectrum, and Qxy(fk)is the quadrature spectrum. 
Subscripts (x,y) indicate this is the cross spectrum between two time series Cj^(t) 

and Cy(t). 

The ordinary coherence function between two stationary records can then be 
determined by 
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and the phase angle by which series x leads series y at the frequency (0 is 
determined by 



Q,y(fl = Tan-i 
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The calculations for the significance level follows Halpem, (1973). A 95 % 
confidence level was used as input into PROSPECT. The significance level (p*) 
represents the coherence value which is exceeded by 5 % of the coherence 
estimates when correlating two random value series. The significance level is 
given by 
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where v is the estimated degrees of freedom and p is the probability such that 
lOO(l-p) % of the coherence values will be greater than p* when the series is 
uncorrelated. 
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V. OBSERVATIONS 



A. SPATIAL AND TEMPORAL CHARACTERISTICS 

To illustrate the vertical structure of the current as it responds to local forcing 
at PI during the deployment of the SC-ADCP, wind stress and current velocities 
are presented in stick vector format with north oriented up (Figure 10). All data 
have been LLP filtered as described in Chapter FV and the local wind stress and 
current from every other bin is plotted as a function of time using six-hourly 
marks. 

The current exhibited a high degree of correlation with the local wind stress. 
Both wind stress and current were strongly polarized in the equatorward 
direction. Equatorward wind stress forcing produced currents whose vertical 
structure was highly sheared with velocities on the order of 35 cm s'^ seen in the 
near surface bins decreasing to 5 - 10 cm s' Mn the bottom bins. There was also a 
pronounced time lag in the response of the water column with depth as the 
effects of wind forcing translated itself deeper into the water column. This was 
readily evident in the 21 - 28 March equatorward wind event. This type of 
response however, appeared to be limited to the first half of the time record. The 
strongest wind event of the record occurred during 24 April - 02 May with wind 
stress magnitudes as great as 3.0 dyn cm'^, however, it did not produce a strong 
equatorward current like the earlier equatorward wind events. Possible 
mechanisms to account for this will be addressed in the discussion section. 
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Figure 10. Time series, in stick vector format, of the low-passed filtered 
local wind stress taken from the Monterey Bay NOAA buoy and current 
velocity from the SC-ADCP at mooring PI: The direction of the vector 
indicates direction of flow, the length of the vector represents magnitude. 
Velocities are in cm s*^, and north is toward the top of the page. 
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Interspersed with the equatorward wind events were relaxations (periods of 
no wind) and wind reversals. There were four distinctive poleward current 
responses in the record which occurred 28 - 29 March, 14-16 April, 02 - 05 May, 
and 08 - 09 May with each corresponding to a relaxation or wind reversal. The 
resulting vertical structure was more uniform than in the equatorward events, and 
did not show any observable time lag in the vertical. This was most clearly seen 
in the latter two events where poleward current magnitudes up to 10 cm s'^ were 
present throughout the column. The response of the water column also appeared 
to be more energetic later in the time series. Extensive research concerning the 
barotropic response of the continental shelf water was conducted during CODE. 
The reader is referred to Kosro (1987), Winant et al. (1987) and Huyer (1987) for 
detailed discussion of similar responses to wind relaxations. 

The V - component plot (Figure 11) clearly showed the vertically sheared 
nature of the equatorward wind stress events and the trend of decreased 
equatorward current with time. The barotropic nature of the poleward relaxation 
events were also evident and as time progressed the tendency for a more 
energetic response was observed. The u - component plot (Figure 12) indicated 
that offshore flow during equatorward events was highly sheared and confined 
to the upper 30 m of the water column while a more constant magnitude (0 - 5 cm 
s'^) onshore "return" flow encompassed the remaining water column to the 
bottom. 

The vertical shear values in cm s*^ were calculated as vertical differences in 
the magnitude of the current at each successively deeper bin (Figure 13). The 
results highlight the shear structure of the v component that resulted from 
equatorward wind forcing and relaxation events. The current response to the 
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Figure 11. Time series of v current component from the SC-ADCP at 80 m 
depth: Velocity values are represented by contours at 5 cm s"^ intervals, as 
indicated by the scale on the right. The contours are estimated from 6 - hourly 
values of 8 - m intervals in the water column. 
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Figure 12. Time series of u current component from the SC-ADCP at 80 m 
depth; Velocity values are represented by contours at 5 cm s"^ intervals, as 
indicated by the scale on the right. The contours are estimated from 6 - hourly 
values of 8 - m intervals in the water column. 



33 






equatorward wind event of 21 - 28 March exhibited shear on the order of 1.5 cm 
s'^ m"^ in the upper depth bins, near zero at mid-depth, and increasing shear to 
the near bottom. In contrast, the shear observed during the relaxation of 02 - 05 
May exhibited a different structure. From the near surface bin, to mid-depth, the 
shear was near zero. It then increased rapidly to a maximum of 1.5 cm s"^ m'^ at 
68 - 72 m, then rapidly to zero. The wind reversal of 14 - 16 April produced a 
current response similar the 02 - 05 May event from mid-depth to the bottom, 
however, a shear of 1.0 cm s"^ m"^ was present in the near surface which was not 
seen in the May relaxation event. 

The u - component shear plot (Figure 14), provides insight on the vertical 
structure of the water column as it responded to upwelling favorable 
equatorward winds. Two equatorward wind stress events that occurred during 
21-28 March and 07-14 April produced large shear, on the order of 1.5 cm s'^ 
m"^ in the near surface, a minimum at mid-depth with shear near zero, and minimal 
shear, less than 0.5 cm s'^ m'^ , from mid-depth to the bottom. 

These shear plots suggested most of the water column was under the 
influence of frictional forces, which induced shear in the water column. These 
forces were acting on both the surface currents, via the surface wind stress, and 
on the near bottom currents, via bottom friction. 

B. EQUATORWARD EVENTS 

To examine the response of the water column to equatorward wind forcing in 
more detail, four daily averaged current profiles were chosen for additional 
analysis (Figure 15). Each profile was computed at least one day into a well 
developed equatorward wind event, and possessed similar characteristics. The v - 
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Figure 13. Time series of vertical shear for v current component: Values of 
shear are calculated by differencing each successive depth bin. The units on the 
y - axis are cm s'^m'^. Vertical lines indicate current events discussed in the text. 
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Figure 14. Time series of vertical shear for u current component: Values of 
shear are calculated by differencing each successive depth bin. The units on the 
y - axis are cm s'^m'1. Vertical lines indicate current events discussed in the text. 
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component in each profile had a magnitude of about 30 cm s'^ at the surface, 
decreasing with depth, and is everywhere equatorward. The temporal trend of 
decreasing equatorward flow was evident when comparing the 7 - 10 March and 
7 -12 April profiles, particularly in the lower portion of the profiles. The current at 
68 - 72 m of 7 - 10 March was on the order of 10 cm s*^ while the magnitude of 
the current on 7 -12 April at the same depth had decreased to near zero. 

Examining the u - component for each event showed that the transition from 
an offshore flow to an onshore flow generally occurred at about 30 m. The 
current was sheared below the transition depth and became less sheared with 
increasing depth. These profiles suggest a wind-driven offshore flow in the 
upper (30 m) layer and a compensatory onshore flow below. Similar results of 
offshore flow in a thin surface layer with onshore flow in a thicker interior layer 
are observed in Huyer (1983). 

The vertical structure of the same four events were projected on a horizontal 
plane as hodographs (Figure 16). The current measured at each bin was 
represented by a dotted vector. The wind stress, represented by a solid vector, 
was increased by a factor of 10 to facilitate comparison with the current. The 
hodographs were similar for all four events. Local wind forcing caused a near 
surface current to develop to the right of the wind stress vector with velocities on 
the order of 20 - 25 cm s' The magnitude of the current was greatest in the near 
surface bins and exhibited a general decrease with increasing depth while veering 
counter-clockwise. The lowest magnitudes on the order of 5 - 10 cm s'^ occurred 
near the bottom. 
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DAILY AVERAGED CURRENT PROFILES 
(EQUATORWARD WIND STRESS EVENTS) 



C COUF (SOUt). r C0MF (OJLSHtS) 




Figure 15. Daily averaged equatorvvard wind forced current profiles In 
component format: The u component of current is represented by a solid line, 

and the v component by a dashed line. Velocity values are in cm s"^ and were 
measured from 14 m to 80 m. 
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Figure 16. Daily averaged equatorward wind forced current profiles in 
hodograph format: The vertical current values at each depth between 14 and 
80 m are projected on a horizontal plane. Current magnitude and direction are 
indicated by dashed vectors. Wind stress is increased by a factor of 10 for clarity 
and is represented by a solid vector. Velocities are in cm s'^. 
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C. POLEWARD EVENTS 

The response of the water column during relaxations and wind reversals was 
significantly different than that observed during equatorward wind forcing. Four 
daily averaged current profiles are presented in the same format as the 
equatorward events, in order to examine the differences in current structure 
(Figure 17). The lack of strong equatorward wind forcing resulted in poleward 
alongshore current with magnitudes ranging from 2 - 3 cm s'^ on 28 - 29 March 
to 10 cm s*^ on 02 - 05 May and 08 - 09 May. The across-shore response 
showed a trend towards increasing offshore flow with increasing depth, which 
was more evident as the strength of the poleward alongshore flow increased with 
time. These profiles suggest the classical case of a uniform upper layer flow 
which backs to the left in an Ekman spiral once the bottom influence is felt. As 
the V - component slowed down due to the effect of the bottom friction, the 
Coriolis force, which is proportional to the current velocity, decreased. The 
pressure gradient which was driving the current is not balanced and the current 
will then swing left (off-shore) until the vector sum of the Coriolis and friction 
forces can balance the pressure induced current (Pond and Pickard, 1983). Using 
the depth of the Ekman spiral as an estimate of the height of the bottom 
boundary, values on the order of 50 - 60 m were inferred. This is consistent with 
the results of Lentz and Trowbridge (1991) who indicate that the bottom mixed 
layer can exceed 50 m in northern California coastal upwelling regions during 
periods of weak stratification. 
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DAILY AVERAGED CURRENT PROFILES 

RELAXATIONS (POLEWARD EVENTS) 
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Figure 17. Daily averaged wind relaxation current profiles in component 
format: The u - component of current is represented by a solid line, and the v - 
component by a dashed line. Velocity values are in cm s'^ and were measured 
from 14 m to 80 m. 



The same events presented in hodograph format (Figure 18), showed no 
relation between the current structure and classical Ekman forcing. Most of the 
flow was oriented poleward, with a decrease in magnitude and backing to the left 
as a function of depth. The surface flow was not oriented to the right of the wind 
stress. The vectors became stronger and more northerly (as opposed to 
alongshore) with time. 
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Figure 18. Daily averaged wind relaxation current profiles in hodograph 
format: The vertical current values at each depth between 14 and 80 m are 
projected on a horizontal plane. Current magnitude and direction are indicated 
by dashed vectors. Wind stress is increased by a factor of 10 for clarity and is 
represented by a solid vector. Velocities are in cm s*T 
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D. COHERENCE AND PHASE 

To quantify the observations of strong visual correlation between alongshore 
wind stress (xy) and the v - component, coherence and phase were calculated for 
Ty and current in both the near surface and the mid-depth bins (Figure 19). 
Strong coherence (above the significance level of 0.73) between Xy and the 12 - 
16 m depth bin were observed at frequencies lower than 0.1 cph. Xy and 
current were in phase at 0.01 cph and as the frequency increased to 0.1 cph, Xy 
leads the current. The coherence between Xy and the v - component at 52 - 56 m 
was similar except for a lack of coherence in the lowest frequency bin. Xy was in 
phase with the current up to 0.07 cph at which point Xy began to lead the 
current. There is significantly less coherence between Xy and the u - component 
at either depth. Coherence at the 12 - 16 m level was only above the significance 
level at 0.9 cph with the current leading Xy. Deeper in the water column the 
coherence between Xy and less coherence between Xy and the u - component at 
either depth. Coherence at the 12 - 16 m level was only above the significance 
level at 0.9 cph with the current leading Xy. Deeper in the water column the 
coherence between Xy and the u - component was high at low frequencies (0.01 
- 0.02 cph), with the u - component roughly 180° out of phase with xy. 

The vertical structure of the current itself was examined with the coherence 
and phase between the shallowest bin at 14 m compared with progressively 
deeper bins. The coherence in the v - component plot (Figure 20) between 14 m 
and 46 m was near unity at all periods longer than about 4 days. Below 46 m 
there was a rapid fall-off in coherence with increasing vertical separation and 
below 62 m the V - component currents were uncorrelated with the 14 m bin. The 
correlation between the 14 m bin and progressively deeper bins in the u - 
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component plot (Figure 21) diminished at much shallower depths and below the 
30 m level, the current was uncorrelated. The results of the spectral analysis 
illustrate that different structure occurred in the upper and lower part of the water 
column, in spite of the lack of density stratification. 
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Figure 19. Coherence and phase between alongshore components of wind 
stress (Ty) and the u and v current components at 12 - 16 m and 52 - 56 m 
depth: The significance level is indicated by a dashed line. 
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is indicated by a solid line. 
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Figure 21. Coherence and phase for u current components at increasing 
depths: The depth bin near the surface (14 m) is compared with each 

successively deeper bin (in 8 m increments down to 78 m). The significance level 
is indicated by a solid line. 
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VI. THEORETICAL MODEL 



To describe the current off Point Sur a theoretical model forced by observed 
wind stress was developed. The model is an adaptation of the elementary wind 
forced current system near a coast model developed by Ekman (1905) and 
Csanady (1982). The following assumptions were made in the model: a long 
straight coast is oriented parallel to the y - axis; steady wind stress blows parallel 
to the coast with no variation in either direction; homogeneous ocean; constant 
eddy coefficient, and water depth greater than the Ekman depth. The 
requirement for a homogeneous fluid seemed justified based on two separate CTD 
surveys conducted on 02 May in the vicinity of the SC-ADCP. These surveys 
showed only weak stratification in the water column (Figure 22). It is assumed 
that surface heating in May was insufficient to stratify the water column in the 
presence of mixing by wind forcing and tidal currents. The assumption of a 
homogeneous ocean then will apply equally well earlier in the record when the 
CTD observations were made. 

Steady equatorward wind stress will deflect the water in the near surface to 
the right, away from the coast. This Ekman transport will create a surface slope 
perpendicular to the coast, and an alongshore flow develops throughout the 
water column which is in geostrophic balance with the across-shore pressure 
gradient force. An adjustment drift in an across-shore u - component allows the 
alongshore flow to develop. As the alongshore current reaches the sea floor, the 
effects of friction will slow the current creating a bottom Ekman layer with 
enhanced across-shore velocities relative to the middle part of the water column. 
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Figure 22. CTD surveys of site PI: CTD casts were made on 02 May in the 
vicinity of PI. 
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Following Ekman (1905), the resulting flow regime will consist of the surface 
and bottom Ekman layers superimposed upon the mostly geostrophic interior 
component (Figure 23a). The characteristics of each flow component can be 
identified in the accompanying hodographs. The wind stress creates a surface 
Ekman regime (labeled a in Figure 23b). The interior component is represented by 
vector (c) and will flow parallel to the coast. The bottom Ekman component (d) 
will develop in response to the frictional drag acting on the interior component. 
When the water depth varies as a function of the Ekman depth (Figure, 23c), 
there is a pronounced effect on the vertical structure of the current. If water 
depth (d) / Ekman depth (D) is less than 0.5 there is a loss of the surface and 
bottom Ekman layers and the total flow in the water column is directed parallel to 
the coast. 

A. GOVERNING EQUATIONS 

Since the interior and bottom components develop because of the surface 
Ekman flow near the coast, the surface component will be derived first followed 
by the interior and bottom components. It is assumed, as in Ekman (1905), that 
the surface layers do not feel the bottom initially, which is the usual frictionless 
development of the problem. It is only some time later, after geostrophic 
adjustment, that bottom friction becomes important. The model approach is 
different in that Ekman (1905) sp»ecified a surface slope rather than allowing the 
wind to develop it. Csanady (1982), however, in allowing the set-up by the wind 
to evolve properly, did not include bottom friction. This model essentially 
combines the two approaches and allows the wind to force the alongshore 
coastal current with the bottom Ekman spiral included as well. 
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Figure 23. Theoretical dynamic coastal current regime accompanying 
alongshore wind: (a) The coastal theoretical wind driven current regime 
showing the u and v - component flow is presented in the top schematic, (b) The 
theoretical current structure is presented (left) in hodograph format, with the 
vectors representing velocity direction and magnitude, (from Neumann and 
Pierson, 1966). (c) A series of hodographs (right) are given where the vectors 
also represent velocity direction and magnitude and where the water depth is a 
function of Ekman depth (from Ekman, 1905). 
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1. Surface Component 

The derivation of the of the surface wind forced component, or surface 
Ekman current, was done by Ekman (1905), and assumes an infinitely deep 
ocean. Ekman (1905), stipulates that the finite depth assumption causes small 
deviations from the mean value of the angle between the wind and the infinite 
depth case, as long as the bottom depth (H) is more than twice the Ekman depth 
(D). The motion occurring can be considered as it would in an infinitely deep 
ocean. Assuming initially a maximum Ekman depth on the order of 30 - 50 m 
(Brink, 1983) and a measured shelf depth of 84 m at mooring PI, the derived 
Ekman current based on an infinitely deep ocean should yield acceptable results. 

As the Ekman equations are well known and widely used in the 
Oceanography community, this section will focus on the assumptions used to 
develop the Ekman equations from the equations of motion following Pond and 
Pickard (1983) and Ekman (1905). Beginning with the equations of motion 



8u 3u 9u 9u Op c 

+ Uz— + V=r- + Wz— = - + fV 

dt dx dy dz dx 



A(a^UA(a^’ 

dxl "dxf dyndy) dz\ "dz, 



(la) 



dv dv 

1— + Uz^ + 

dt dx 



dv dv dp ^ 

^ + w— = - - fu 

dy dz dy 



+ 




(lb) 
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where Ax, Ay, and Az are the respective eddy coefficients with the specific 
volume anomaly a = 1 / p. 

Equations (la), (lb), and (Ic) contain non-linear terms and generally are 
not explicitly solvable. By scaling these equations utilizing rough estimates of 
the size of each, neglecting the non-linear terms, and only considering frictional 
forcing acting in the vertical direction, it is possible to neglect components of 
equations (la), (lb), and (Ic) leaving the remaining terms to describe the actual 
motion. Considering only the horizontal equations of motion, equations (la) and 
(lb) become. 



f. A 

fv = a5^-A2 — 
ox dz2 



( 2 ) 



fu = - Az 

ax az2 



(3) 



There are two causative forces for motion in equations (2) and (3), the 
distribution of mass (density) which gives rise to pressure terms and the wind 
friction terms. If the velocity is considered having two parts, one associated with 
the horizontal pressure gradient and the other with momentum transferred by 
vertical friction, each can be solved for separately and then reconstructed. 
Taking the u - component equation (2) 
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where 



and 



fvg = 

ox ug, vg are Ekman velocity components 



fve = - Az 



9^(Ug + Ue) 
0z2 



Ue, ve are Ekman velocity components 
associated with vertical shear faction 



As the water is assumed homogeneous and there is no initial surface 
slope, the pressure terms are equivalent to zero and therefore ug and vg are zero. 
Similar arguments apply to the v - component equation. The resulting Ekman 
equations are: 



fve+ Az^-^ = 0 
3z2 


(5) 


-fue+ Az^-^ = 0 
3z2 


(6) 



The solutions to equations (5) and (6) are well known (Ekman, 1905; 
Neumann and Pierson, 1966; Pond and Pickard, 1983) and yield 



(7) 



where 



u = Voe-(S)^cos(i-gz) 

v = Vt,e'(D)^s'n(5--0z) 



D = Ekman 



Depth = 71 



V 



MI 

pf 



( 8 ) 

(9) 



and 

Vo = Surface Current Velocity = —J=z= 

^ pA^f 

Equations (7) and (8) are the u and v - component currents for the 
surface Ekman layer. From these equations it is observed that at z = 0, 

u = Vo cos 45° 

V = Vo sin 45° 

or that the surface current will develop oriented 45° to the right of the stress in 
the Northern Hemisphere. By varying z the Ekman Spiral will be produced for 
the pure wind driven current. 

The Ekman current has a maximum speed at the surface and decreases 
with increasing depth. If equations (7) and (8) are integrated over the Ekman 
depth (D), the important transport equations are given as 

and Ty =0 (jo) 
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All the transport in the Ekman layer due to the alongshore wind stress is in the 
across-shore direction. This has important implications for coastal upwelling as 
described next. 

2. Interior Component 

Surface Ekman drift in the presence of a coast will create pressure forces 
in the interior of the fluid. Following the derivation of Csanady (1982), the 
interior component will be shown to consist of an alongshore geostrophic flow 
and an "adjustment drift" current in the across-shore direction. The horizontal 
transport equations are 





(11) 




(12) 


au av ac 
ax at ' at 


(13) 



where c^ = gH, (F^, Fy) are the kinematic surface wind stress components, and 
(Bjj, By) are the kinematic bottom stress com.ponents. 

Differentiating (1 1) with respect to y, (12) with respect to x, subtracting 
(12) from (11), and substituting (13) into the resulting equation yields 



_a 

at 




(14) 
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This is the depth integrated analog of the linearized potential vorticity equation. 
With the spatial variability of the wind stress and bottom stress ignored, and the 
motion starting from rest, (z = 0 at t =0), equation (14) becomes 

dx * dy (15) 

Differentiating (11) with respect to x and (12) with respect to y, summing and 
substituting into (13) yields 



^-c2 V^^ + fX = 0 
dt^ 



(16) 



Transports are now found by assuming a coast parallel to the y - axis, 
with alongshore wind in the (+) y direction imposed at t = 0. Neglecting bottom 
friction, the boundary condition at t = 0, is 

Fx = 0 Bx = 0 
Fy = u*^ By = 0 



Substituting into equations (11) and (12) 



dt ox 



(17) 



dt 



+ fU = U*2 



(18) 
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Imposing the boundary condition at the coast, U = 0, and starting the motion from 
rest equation (17) and (18) yield 



V = U*2 t (x = 0) 


(19) 


dx c2 


(20) 



Equation (19) represents the flow that would result if no rotation were present. 
The fluid near the coast simply accelerates downwind. Equation (20) shows that 



the developing sea surface slope ^ is a function of the wind duration, magnitude 

dx 

and the free wave speed c^ = gH. 

To establish the pressure field, the full solution of the impulsively forced 
problem, equation (16), is found using Laplace Transform methods (Crepon, 1967) 
the result for the free surface elevation is 



C = u 





( 21 ) 



where J© is a Bessel function and Y(X) is the Heavyside function (Y = 1, X > 0, Y = 
0, ^ < 0). A long time after the imposition of the wind stress equation (21) may be 
expressed as 
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y _ U*^ 

^ fc 



(ft)ex/R-^Icosl-^ 
n vFt 



ft»l 



( 22 ) 



steady 



oscillatory 



which is the pressure field described by a flow confined to a coastal band of scale 
width R. The steady amplitude will increases linearly with time and the 
oscillatory part will decay with time. 

The non-oscillatoiy' part of (22) is; 

u^2 t gx/R 

^ < 23 ) 



where x is the distance from the coast and R is the external Rossby radius of 
deformation given as 



R = 



VgH 

f 



The alongshore transport is found by differentiating ^ with respect to x 



3x R 




(24) 



and substituting into equation (17) at x = 0 and U = 0. 

This yields V = u*^ t e’^ (25) 

where R = c/f. Using equation (25) and (18) gives the integrated across-shore 
transport 
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(26) 



U=M^(1 - e^) 



Equations (23) and (25) with 



fV = c 



2 ^ 

dx 



show that the alongshore transport V is in geostrophic balance with the pressure 
field, not only at the coast but everywhere in the coastal band of scale width R. 
Far from the coast, V goes to zero as the sea surface goes to zero. There is also 
significant across-shore transport U, over most of the coastal band except very 
near the coast where x/R «1. The U transport builds up according to the time 
scale f . which at this latitude is 20 hours, and since there is no along-shore 
pressure gradient, U is not in geostrophic balance. Near the coast where U is 
approximately zero the wind will only accelerate the water alongshore. Far from 



the coast U = 




the Coriolis force associated with the Ekman drift balances 



the wind stress and only surface, Ekman flow will develop. In between, the flow 
transitions smoothly between the two states (Figure 24). 
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Figure 24. Coastal transport schematic: Alongshore transport (V) is in 

geostrophic balance within the coastal band from the coast to - x/R » 1. There is 
also significant across-shore transport (U) over most of the coastal band except 
very near the coast where -x/R « 1 , 
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The solutions to the previous equations have yielded the local sea level 



gradients 



^ K 

ax ’ dy 



which are the depth independent forcing terms of equations (16) and (17). 

At this point the interior velocity distribution will be considered and the 



local solutions examined for equations (16) and (17), 
v(z,t) which depend parametrically on x are written as 


These solutions u(z,t) and 


u = ui(t) -1- U 2 (z,t) 


(27) 


V = Vi(t) + V 2 (z,t) 


(28) 


where ui(t)and vi(t) is a solution to equation (16) and (17) with the stress 


gradient terms deleted. 

3ui n 

at -f''' = -8ax 


(29) 


0Vi r 

a. ^f“' = -«ay 





Earlier it was shown that the global solution yielded depth-integrated 
transports U and V. Equations (27) and (28) will also apply to the transports, so 
with the components ui(t)and vi(t) independent of depth, the equations for the 
interior velocity distribution become 



(32) 



V]=jL(v.V2) 



To solve for the interior velocity, two steps are involved. The 
frictionally induced transport (U 2 and V 2 ) must be calculated. The pressure 

induced interior velocities are then solved by subtracting (U 2 and V 2 ) from the 

total transport components (U and V). The frictionally induced transport 
components (the Ekman transport), are derived from the global solutions to 
equations (16) and (17) as 



fV = 0 
fU = F 



3x dy 



= 0 



which are satisfied by 

U2- — 

V2 = 0 

where F = 

Solving for Uj and Vj from equations (29) and (30) yield the pressure 
field induced velocities 
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( 33 ) 




;X/R 



Vi = e’^ 

H 



(34) 



Equations (33) and (34) are functions of exp(x/R) which determines the coastal 
decay scale length. This is the region where wind stress acting alongshore 
indirectly produces a current in its own direction. For site PI, moored 5 km off 
the coast, exp(x/R) = exp(5/330) = 1.0, so PI is well within the coastal band 
where these dynamics apply. Therefore Uj is simply the adjustment drift and Vj is 

the geostrophic transport generated by the wind stress acting over time (t). 

dC 

Referring to equation (30) and taking r^ = 0, the alongshore geostrophic 

dy 

component must be generated by an appropriate across-shore motion during a 
spin-up period which occurs on a time scale f This equates to a spin-up time 
on the order of about 20 hours at 36® 20 ' N latitude. 

After initial set-up, the magnitude of Vj will increase as long as forcing 

continues. In the frictionless case, the current acceleration can be infinite 
however, when bottom friction is introduced, the surface stress will at some point 
be balanced by bottom stress and the current will no longer accelerate. In this 
case the bottom stress (By) will be equivalent to surface stress (xy) and is given 
by 



Xy = By = V^ Cd 



(35) 
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where v is the near bottom velocity in the alongshore direction and Cd is the drag 
coefficient. Observed LLP values of the v - component in the deepest depth bin 
are on the order of 10 cm s'^ , however, to obtain a more accurate estimation of 
the V - component, including tidal contributions, unfiltered (15 min) current data 
were analyzed. Typical near-bottom maximum v - component values of order 20 
cm s'^ were observed throughout the record. This range of velocities, when used 
with a surface stress of 1 .0 dyn cm'^ in equation (35) yields a range of values for 
Cd from 2.5 - 10 x 10'^. These values are slightly higher than the Cd used by 
Trowbridge and Lentz (1991), but within the range of values reported in Table 2 
of Grant et al. (1983). 

3. Bottom Component 

In order to develop a more realistic model of the shelf flow, a bottom 
Ekman spiral is allowed to form (Ekman, 1905; Pond and Pickard, 1983) 
Assuming a constant Az, the current flowing over the sea floor will generate an 
Ekman spiral above the bottom to a height D equal to the depth of the surface 
Ekman spiral however, the bottom spiral backs counter-clockwise with depth, 
rather than clockwise. The Ekman equations (5) and (6) will still apply but the 
boundary conditions will be different. At the bottom the velocity must vanish, 
(u = 0, V = 0) and must reach a constant value some distance above the frictional 
layer. Assuming this layer, above the frictional layer, is in geostrophic balance it is 
then independent of z. Taking u = 0 and v = vg in the geostrophic region the 
solution to (5) and (6) with the new boundary currents gives 

Ub = [{e-"z/D) sin (tiz/D)] 

H 
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Vb = ^ cos {tizA>)] 

H 



where D is the bottom Ekman layer height defined by (9) as before and z is now 
the height above the bottom. 

4. Solutions 

The total solution will be the sum of the surface, interior, and bottom 

layers: 
u(x, z, t) = 

Voe-(S)"cos(S--iz) - e*/R ■ z)®) yn (^(H - zjofl , _ 

'4 D / H H (3oa) 

(surface Ekman) (adjustment drift) (bottom Ekman) 

v(x, z, t) = 

Voel5)"sin(5-gz) - e-'/R - ^Vd) cos (7t(H - 

(surface Ekman) (geostrophic) (bottom Ekman) 

The u and v - component equations will consist of contributions from both 
surface and bottom Ekman components. The u - component will also contain an 
adjustment drift contribution and the v - component, a geostrophic contribution 
in the interior. Where the bottom Ekman part has been re-written so that z is 
positive downward from the surface as it is in the surface Ekman layer. 
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B. RESULTS 

The initial input variables for the governing equations were chosen based on 
observations taken from the 08 - 09 March equatorward wind event. The surface 
stress is computed by averaging the wind stress between 08 and 09 March. The 
initial choice of Az is computed from the solution of (9) using a value of D 
estimated as the zero crossing point (about 30 m) in the u (across-shore) velocity 
component. This is based on two-dimensional upwelling theory which stipulates 
that offshore transport due to alongshore wind stress must be balanced by 
onshore transport below the Ekman layer (Brink, 1983). Using this procedure a 
value of Az = 4.5 x 10'^ was computed. The initial duration of the wind event 
was chosen as 20 hours (the required initial spin-up time) but longer values may 
be appropriate. 

The hodograph of the resulting theoretical vertical velocity field clearly 
shows the current structure of the three flow regimes (Figure 25). With the wind 
stress acting along the (-) y axis parallel to the coast, a near surface Ekman layer 
develops to the right of the wind stress vector. As this Ekman layer is 
superimposed on the alongshore pressure driven velocity field, the magnitude of 
the current does not go to zero with depth rather, it approaches the magnitude of 
the geostrophic interior current. The usual clockwise turning of the Ekman spiral 
is also suppressed by the geostrophic component. Thus the uniform interior 
current is represented by a bunching of the current vectors at mid-depth in the 
general direction of the wind stress at a magnitude of 20 - 25 cm s*^. As the 
interior current begins to feel the bottom, friction causes a decrease in magnitude 
and the current vectors respond to this and veer counter-clockwise. The deepest 
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vector plotted is at 78 m to match observations, although the theoretical 
hodograph does back to zero at the bottom. 

The theoretical results are compared with the observed results in the 
component plot (Figure 26). The model reproduced the general structure of the 
current quite well. The v - component is equatorward throughout the water 
column and the u component transitions from off-shore to on-shore in the near 




Figure 25. Theoretical current structure: Current direction and magnitude, 
at depth increments of 8 m starting at the surface, are represented by dashed 
vectors and presented in hodograph format. The wind stress is indicated by a 
thick vector pointing down with the magnitude increased by a factor of 10 for 
comparison. The units of wind stress are dyn cm*^, and the units of velocity are in 
cm s"^ . 
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Figure 26. Theoretical current vs observed current: The theoretical current 
is based on inputs taken from the 9 March observations and is compared with the 
observed current structure on the same date. The u - component is represented 
by a solid line the v - component by a dashed line. Theoretical current 
components are plotted from the surface to 80 m (8 m increments). 
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surface, both with about the right magnitude. The u and v - components of the 
model however, do not possess all the characteristics of the observed flow. The 
observed u and v - components do not possess the constant magnitude interior 
regime that is present in the model, and both are sheared in the interior portion of 
the water column. 

The initial solution of the model was based on an estimate of a constant Az. 
As Az is a property of the ambient conditions of the water column rather than the 
fluid itself, the initial estimate may vary. Complicated bathymetry in the vicinity 
of PI and the lack of stratification in the water column indicate that a larger value 
of Az may be more appropriate. Utilizing Az as a tuning parameter for the model, 
the solutions to the governing equations (36a and 36b) were repeated for 
different values of Az, all else being the same. Figure 27 is a component plot of 
the theoretical current showing the effect of increasing Az from 2.1 x 10*^ to 26.4 
X 10"^ m^ s‘^. By increasing Az, more velocity is diffused into the interior with an 
associated loss of the constant velocity interior signatures of the u and v 
components observed in the initial model. Note also that the transition depth 
between off-shore and on-shore flow stabilizes at 25 - 30 m, which means that 
this cross over point is actually a poor indicator of the Ekman depth (D) when 
bottom stress is included in the model and the water is not very deep. 

A value of Az of 12.6 x 10 m^ s*^ was chosen as input for the best-fit 
solution because it most accurately represented the u - component and greater 
values of Az caused shoaling of the transition level between off-shore and on- 
shore flow. Substituting this value of Az into the governing equations results in 
the hodograph presented in Figure 28. The vectors now back smoothly counter- 
clockwise with depth as the surface and bottom Ekman layers merge into one 
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Figure 28. Theoretical current structure: The theoretical current is solved 
for based on inputs taken from 9 March observations utilizing a value for Az = 
12.6 X 10-3 m2 s-1 . 
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another. The current structure resembles Ekman's theoretical hodograph for a 
water depth = 1.25 D, (Figure 23c) The plots of the theoretical and observed 
current for 9 March (Figure 29) indicate that the larger value of Az provides a 
better approximation of the actual current. The transition depth and general 
shape of the u - component in the theoretical current is approximately the same as 
the observed current. The vertically integrated across-shore transport is not 
balanced, however, it is not used as a criteria for model comparison because: 1) it 
is quite sensitive to the choice of coordinate system, and 2) it is not clear that the 
observed across-shore transport should necessarily be balanced near a coastal 
promontory such as Point Sur. Cold upwelling filaments are often observed there 
(Breaker and Mooers, 1987) which presumably requires alongshore convergence 
or divergence and an associated across-shore transport. The theoretical v - 
component more closely resembles the profile of the observed current than the 
curve found using Az = 4.5 x 10'^ , but still exhibits less shear than the observed 
current. 

Utilizing Az = 12.6 x 10 m^ s‘^ and taking the duration and magnitude of 
the wind stress from each of the four equatorward events as inputs to the 
governing equations, the other equatorward wind stress events were calculated 
(Figure 30). The general characteristics of each event are reproduced well. The 
nature of the theoretical plots is such that velocity components up to the surface 
are calculated. As a result the model not only depicts the current structure but 
also has a hindcast capability for filling in the missing data in the near surface 
bins. The accompanying hodographs of the theoretical and observed currents 
(Figure 31) also illustrate this capability. The model predicts surface values 
ranging from 20 - 35 cm s'^ for these wind events. 
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Figure 29. Theoretical current vs observed current: The theoretical current 
is based on inputs taken from the 9 March observations with an increased value 
for Az = 12.6 X 10'^ m^ s"^ substituted.and is compared with the observed current 
structure on the same date. The u component is represented by a solid line the v 
component by a dashed line. Theoretical current components are plotted from 
the surface to 80 m in 8 m increments. 
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Figure 30. Theoretical current vs observed current: Theoretical current is 
based on inputs taken from the equatorward wind stress observations and 
compared with the observed current values from the same time frame. The u 
component is represented by a solid line the v component by a dashed line. 
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Figure 31. Theoretical current vs observed current hodographs: 
Theoretical current is based on inputs taken from the equatorward wind stress 
observations and compared with the observed current values in a series of 
hodographs from the same time frame. 
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VII. DISCUSSION 



A. MODEL RESULTS 

1. Features 

The simple wind driven coastal model reproduced the main features of 
the coastal current regime surprisingly well considering the complexity of the 
measurement site. It was particularly adept at describing the cross-shelf 
component of flow. The modeled transition from offshore to onshore flow 
occurred at approximately the 30 m level and the magnitude of the adjustment 
drift, at 10 cm s' ^ , was very similar to the observations. The magnitude of the 
interior alongshore component was sensitive to the duration of the wind stress. 
Using t = 24 h in the model produced surface (20 - 25 cm s'^), and interior 
velocities that agreed with the observations. Beyond this, the surface stress was 
apparently balanced by the bottom stress and further accelerations did not occur. 
The interior portion of the v - component of the model was constant for all the 
equatorward wind events. Also, the constant magnitude interior flow was more 
pronounced when using smaller values of the wind stress. The sheared interior 
flow in the v - component could not be reproduced without using Az coefficients 
that were too large to be realistic. 

2. Causative factors for observed model discrepancies 

The discrepancies existing between model results and the observed 
current may possibly be due to variations in the assumptions utilized in deriving 
the governing equations of the model. Physical factors which were not included, 
but which may be important are discussed. 
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a. Slight baroclinic effects and variable Az 

These two effects are actually interelated. Slight stratification 
causes no appreciable difference from the homogeneous case (Neumann and 
Pierson, 1966). A well mixed upper layer over a strong thermocline responds 
differently, primarily due to a much greater Az in the well mixed layer than in the 
stratified and lower layer below. In this case, the upper layer will "slide" over the 
lower layer, and the surface velocity and angle of deflection will be greater than 
in the unstratified case. No evidence was seen for this in early spring off Point 
Sur, and the assumption of a homogeneous ocean is believed to be good. 

b. Nonuniform wind stress 

The presence of spatial variability in alongshore forcing (xy) imply 
the assumptions made in equation (9.5) 

Fx = 0 

Fy = U*2 

and utilized in developing the interior current components along an infinitely 
long straight coast were not strictly valid. Alongshore variation in Xy will caused 

alongshore pressure gradients to develop and introduce varying v - component 
forcing (Fy u*^). These variations influence the observed current structure and 

cannot be reproduced by the simplified model. The effects of a non uniform-wind 
stress field will be discussed in detail in section B. 

c. Bathymetric considerations 

The continental shelf off Point Sur is neither straight or flat. The 
bottom depth shoaled to 55 m depth just 0.5 km to the southeast, in a rough- 
bottomed region known as "the rock pile." The bottom contours were not 
parallel to the coastline, making a choice of coordinate systems difficult. The 
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coordinate system was selected with u - east and v - north, recognizing that true 
"alongshore" flows, in the large-scale sense, would be partially mapped into the u 
- component. This affected the u and v - component plots but not the magnitude 
or sense of rotation of the hodographs. Across-shore transports were very 
sensitive to this problem and are not calculated. Considering the simplistic case 
of continental shelf bathymetry increasing with depth as a function of distance 
from the coast, the across-shore adjustment drift would also experience changes 
in vorticity not accounted for in the model. 
d. Remote influences 

The contributions of remote influences on the observed currents in 
the vicinity of Point Sur were not included in the model results. These influences 
include the presence of long, wind forced, coastal trapped waves propagating 
poleward along the shelf and slope; mesoscale variability over the shelf; and the 
presence of alongshore pressure gradients due to factors other than coastal 
trapped waves. A unique feature of these data relative to the CODE results is 
that most of the current variability can be accounted for via the local, rather than 
remote wind forcing. This is perhaps because the short record obtained from PI 
lends itself to analysis of the higher wind-forced frequencies, which tend to be 
more locally forced (Denbo and Allen, 1987; Winant et al., 1987). There may also 
be a significant difference in the alongshore scales of the wind field between 
CODE and Point Sur regions. To be apparent, remote forcing must necessarily be 
occurring more than a coherence scale away from the observation site (Davis and 
Bogden, 1989). Off Point Sur, this coherence break in the wind field occurs at 
Point Conception, a sharp change in coastline orientation, which may also 
represent a dynamical break in the alongshore current structure as well. 
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Most of the mesoscale variability in the CCS is topographically 
constrained to the continental slope region and farther offshore, and does not 
impact the mid-shelf region. While there are exceptions (Lentz, 1987), mesoscale 
eddies were apparently not important during the short deployment of the SC- 
ADCP off Point Sur. 

Other causes of alongshore pressure gradients include variability in 
the alongshore component of the wind stress, differential upper layer heating, and 
large-scale oceanic influences. These factors are discussed subsequently in 
connection with an anomalous equatorward wind event which did not elicit the 
expected response in the SC-ADCP currents. 

B. ANOMALOUS EQUATORWARD WIND EVENTS 

To highlight the difficulty in modeling the current structure on the 
continental shelf using a simple wind forced model, this section focuses on the 
unusual response to an equatorward wind event during 24 April - 02 May. The 
magnitude of this event was about 3.0 dyn cm"^, however, the current response 
was different from similar wind forced events observed elsewhere in the data set. 
To examine the respective current structures, enlarged stick vector plots for this 
event and an earlier equatorward wind event during 20 - 26 March are shown 
(Figure 32). The 20 - 26 March event has a large near surface magnitude (35 cm 



80 



E 

0 
\ 

1 

iS . 

q: 



o 

z 

Si 



20 - 
0.0 - 
• 2.0 - 



J I I I I I ^ L 



I \\\\\\\\\\\\\\W^^ I 



"H — I — I — I — I — rn — T“ 
19 20 21 22 23 24 25 26 27 

MAR 




I 

20. 7 
0 

-20. I 



20 7 
0 <N 

-20 t 



20 . ? 



-20. i 



E 

20 2 



- - 20 . 



E 

20 « 
h 0. S 



MAR 



CN 





APR MAY 



Figure 32. Time series showing equatorward wind stress events during 20 - 
26 March and 24 April - 02 May: The vector direction and length correspond 
to current direction and magnitude respectively with north up. Velocities are in 
cm s‘^ . 
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s'^X and a highly sheared response, with the current magnitude decreasing to 5 
cm s'^ at 76 - 80 m. The vertical structure is time-lagged with the surface currents 
leading the bottom current. In contrast, surface currents during the 24 April - 02 
May event were 15 cm s' ^ and did not j)enetrate beyond 60 - 64 m. Below this, 
currents were weak (10 cm s'^) and variable, often near zero. 

L Poleward pressure gradient 

A poleward pressure gradient growing stronger with time is a possible 
mechanism to account for the lack of current response to the large equatorward 
wind event of 24 April - 02 May. To examine this hypothesis the difference in 
SSP between Monterey Bay and Port San Luis (ASSP) is presented along with 
the local wind stress data from the Monterey Bay buoy in Figure 33. A linear 
least squares fit to the ASSP data is indicated by a dashed line. A trend of 
increasing ASSP is observed over time. Of particular note is the large magnitude 
(-I-) ASSP that occurs in conjunction with the equatorward wind event of 24 April 
- 02 May. The positive ASSP indicates higher sea level at Port San Luis than at 
Monterey Bay. This infers a poleward pressure gradient which will oppose the 
equatorward alongshore wind stress. 

2. Alongshore variability in the alongshore wind stress 

The observed increase in the poleward pressure gradient may be 
attributed to 1) alongshore variability in Xy, 2) greater upper layer heating to the 

south than the north as the spring season arrives, 3) the presence of coastal 
trapped waves at various frequencies. The anomalously large wind event of 24 
April - 02 May suggest that this particular event may be related to a 
geographically localized synoptic scale event. If so, a large gradient in along- 
shore wind stress could be generated and subsequently contribute to the set-up 
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ADJUSTED SEA LEVEL HEIGHT = PSL-MRY 




Figure 33. A SSP vs wind stress: A SSP is computed as the difference in SSP 
between Port San Luis (PSL) and Monterey Bay. Positive values indicate higher 
SSP at PSL. Differences are measured in cm. 
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of a poleward pressure gradient. This would violate one of the fundamental 
assumptions of the model, of constant alongshore wind stress, and would prohibit 
the model from reproducing the results. 

Alongshore wind stress data from the Monterey Bay NOAA buoy and 

the Santa Maria NOAA buoy, located 200 km to the south near Port San Luis are 
presented in Figure 34. Taking the difference in Xy (Santa Maria - Monterey Bay) 

and superimposing ASSP (Figure 35) shows that the greatest magnitude 
difference in Xy occurs in conjunction with the largest positive ASSP, which 

supports the idea that the 24 April - 02 May wind event had a smaller alongshore 
extent than the earlier wind events. Other variations in ASSP, including the 

3Xy 

increasing trend are unrelated to and must be due to other processes. 
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Figure 34. Alongshore component of the wind stress at Monterey Bay and 
the wind stress at Santa Maria: The wind stress is measured in dyn cm'^. 
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Figure 35. dxyfdy vs A SSP: dxy/dy is the difference between the alongshore 
component of the wind stress at Monterey Bay and at Santa Maria. Positive 
values indicate stronger wind stress at Monterey Bay. 
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VIII. 



CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The purpose of this study was to observe the circulation on the continental 
shelf off Point Sur, California using a SC-ADCP. The SC-ADCP proved to be a 
reliable and accurate instrument for measuring current data. The instrument 
retrieved 100 percent of the available current data from deployment on 28 
February through recovery on 11 May 1990. 

The SC-ADCP was moored on the 84 m isobath with the transducer head at 
80 m. The mooring was on a relatively flat area, h cwever, a substantial rock pile 
shoaling to 55 m was 500 m to the southeast creating very complicated 
bathymetry. An ensemble of 170 one second pings were averaged every 15 
minutes and recorded on a cassette tape. At a vertical resolution of 4 m, twenty 
depth bins were computed per ensemble. Prior to averaging, each ping was 
converted into earth coordinates thus producing three velocity components. The 
mooring geometry provided a stable platform for the SC-ADCP. As a result of 
negligible tilt, current data accuracies of 0.5 - 1.0 cm s'^ were attained to the near 
surface. Due to side lobe bias introduced at the near surface, the top 15 m of the 
current profile was discarded. 

Low passed filtered data indicate high coherence between local wind forcing, 
taken from the Monterey Bay NOAA buoy, and current measured at the SC- 
ADCP. During equatorward wind events, current response was vertically sheared 
and lagged with near surface magnitudes up to 35 cm s' A notable exception 
occurred late in the record when the largest magnitude wind event (3.0 dyn 
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cm ) failed to produce a comparable response in the water column. Analysis of 
the V - current component indicated a trend toward decreasing equatorward flow 
throughout the record. Adjusted sea level (ASL) differences suggested the trend 
may be due to a building poleward pressure gradient as evidenced by a higher 
ASL south of PI. Relaxations and wind reversals produced more barotropic 
poleward responses on the order of 10 cm s*^ in the water column. Possibly due 
to an increasing poleward pressure gradient, these responses are observed to 
increase in frequency and magnitude as the record progresses. 

A simple wind forced theoretical model adapted from Ekman (1905) and 
Csanady (1982) was used to describe the flow in the coastal regime. The model 
.onsisted of surface and bottom boundary layers superimposed on a mostly 
geostrophic interior flow driven by set-up and set-down of the sea surface near 
the coast. Essential features of the model including the magnitude, direction, and 
sense of rotation of the current vectors with depth were reproduced well when 
forced by the actual surface wind stress and realistic vertical eddy coefficients. 

Differences existing between model results and the observed current may 
possibly be due to variations in the assumptions utilized in deriving the 
governing equations of the model. Remote influences can also affect and 
influence the observed currents and thus not be reflected in the model results. 

B. RECOMMENDATIONS 

This study utilized SC-ADCP data to describe the current structure of the 
flow along the central California continental shelf. The SC-ADCP at PI provided 
exceptional temporal and vertically spatial coverage, however, at the expense of 
horizontal coverage. A better understanding of the horizontal spatial variability 
of the shelf flow could be obtained by mooring several SC-ADCPs on the 
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continental shelf or scheduling concurrent hydrographic and shipboard surveys 
in the vicinity of the mooring. It would also be useful to moor the instruments for 
more than four months, to compare the statistical confidence in the wind-forced 
events and allow the discrimination of lower frequencies. This would also allow a 
more rigorous analysis of remote forcing from the south. Additional 
measurements in other seasons would allow analysis of summer stratification, the 
fall transition, and the onset of winter cooling. 

Regarding the positioning of the SC-ADCP for future deployments, 
consideration should be given to the alignment of the instrument along constant 
lines of latitude as done with established moorings P2 and P3. As the vertical 
structure of typical continental shelf currents is likely to change most rapidly in 
the across-shelf rather than the along-shelf direction (Beardsley and Lentz, 1987) 
the alignment of the instruments along the same latitude would aid in the analysis 
of the across-shelf variability in flow structure. 
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